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^ c| We discuss the most important results which are obtained from a 
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recent calculation of the order corrections to the coefficient func- 
tions contributing to the fragmentation functions Fk{x, Q 2 ) (k = T,L,A) 



in e + e -annihilation. Prom these fragmentation functions one can de- 
ed 1 rive the corresponding integrated transverse (cft), longitudinal (cfl) 

and asymmetric (a a) cross sections. The sum a tot = ct + ""L cor- 
rected up to order a 2 agrees with the well known result in the literature 
providing us with an independent check on our calculation. It turns 
out that the order a 2 corrections to the transverse and asymmetric 
quantities are small. However these corrections to Fl{x,Q 2 ) and ctl 
are large so that one gets a better agreement between the theoretical 
predictions and the data obtained from the LEP experiments. 
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The measurement of the fragmentation functions in the process 



e + + e" ^ 1 ,Z^H+ "X", (1) 

provides us, in addition to other experiments like deep inelastic lepton-hadron 
scattering, with a new test of scaling violations as predicted by perturba- 
tive quantum chromodynamics (QCD). Here "X" denotes any inclusive final 
hadronic state and H represents either a specific charged outgoing hadron or 
a sum over all charged hadron species. This process has been studied over 
a wide range of energies of many different e + e~-colliders. The most recent 
data are coming from LEP (see , @] , ) . The unpolarized differential cross 
section of the above process is given by 

d 2 H 3,. 2n\ d °T 3 • 2n d(7 L 3 n d(J A ^ 

— = - 1 + cos 2 6)-/- + - sin 2 + - cos0-A. (2) 

dx a cos v 8 dx 4 dx A dx 

The Bj0rken scaling variable x is defined by 

x = tI' q 2 = Q 2 >0, 0<X<1, (3) 

Q 

where p and q are the four-momenta of the produced particle H and the 
virtual vector boson (7, Z) respectively. In the centre of mass (CM) frame of 
the electron-positron pair the variable x can be interpreted as a fraction of 
the beam energy carried away by the hadron H . The variable 9 denotes the 
angle of emission of particle H with respect to the electron beam direction in 
the CM frame. The transverse, longitudinal and asymmetric cross sections 
in (|) are defined by , , and respectively. The latter only shows 
up if the intermediate vector boson is given by the Z-boson and is absent in 
purely electromagnetic annihilation. 

Before the advent of LEP1 the CM energies were so low (y/s < Mz) that a a 
could not be measured and no effort was made to separate a? from so that 
only data for the sum of the above cross sections in eq. @j were available. Re- 
cently after LEP1 came into operation ALEPH and OPAL obtained data for 
ot and separately and the latter collaboration even made a measurement 
of a a for the first time. The separation of ot and is important because 
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the latter cross section enables us to extract the strong coupling constant 
ot s and allows us to determine the gluon fragmentation density D^(x,fi 2 )) 
with a much higher accuracy as could be done before. Furthermore the mea- 
surement of oa provides us with information on hadronization effects since 
the QCD corrections are very small. In the QCD improved parton model, 
which describes the production of the parton denoted by I (I = q, q, g) and its 
subsequent fragmentation into hadron H, the cross sections oj? (k = T, L, A) 
can be expressed as follows 



dx 



I ^ [^l(Q 2 ) \d$ (~ /i 2 ) Cjg>, Q 2 /» 2 ) + Df g, ^ 



for k = T,L. In the case of the asymmetric cross section we have 



dx 



1 dz 

r Z 



(4) 



(5) 



In the formulae (^) and ([|) we have introduced the following notations. The 
function D^z,^ 2 ) denotes the parton fragmentation density corresponding 
to the hadron of species H. Further we have defined the singlet (S) and non- 
singlet (NS, A) combinations of the quark fragmentation densities. They are 
given by 



D H s {z^ 2 ) = -EK(^ 2 ) + ^(^ 2 )), 

U f q=l 



D"(z^ 2 ) + D*(z,v 2 )-D*(z^ 2 ) } 



D H A (z^ 2 ) = D«{z^ 2 )-D«{z^ 2 



(6) 

(7) 
(8) 



The index q stands for the quark species and nf denotes the number of light 
flavours. Further the variable fi appearing in Df 1 ^,/! 2 ) and the coefficient 
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functions Ck,i(z, Q 2 / (i 2 ) {k = T,L,A; I = q, g) stands for the mass factoriza- 
tion scale which for convenience has been put equal to the renormalization 
scale. The pointlike cross section and the asymmetry factor of the 
process e + + e~ — > q + q can be found in eq. (10) and eq. (15) of [|J 
respectively. The total cross section, summed over all flavours is given by 
cx t ( °|(Q 2 ) = J2gLi ^ 0) (Q 2 )- The QCD corrections in eqs. (§, (§) are described 
by the coefficient functions C r k l (k = T, L, A; I = q, g) which can be distin- 
guished with respect to the flavour group SU(nf) into a singlet (r = S) and 
a non-singlet part (r = NS). The gluonic coefficient function only receives 
contributions from flavour singlet channel partonic subprocesses so that we 
can drop the superscript S on Ck, g . However the quark coefficient functions 
get flavour singlet as well as flavour non-singlet contributions. Up to first or- 
der in the strong coupling constant a s it turns out that = Cf q . However 
in higher order both quantities start to deviate from each other. Hence we 
define the purely singlet coefficient function C£g via 

Ck,q = Ck,q + ^k,q • (9) 

Like Ck, g the purely singlet coefficient function only receives contributions 
from the flavour singlet channel partonic subprocesses which in the former 
case show up for the first time in order a 2 . 

Using charge conjugation invariance of the strong interactions one can show 
that Cj[ t q = —Ca^ an d Cji S q = ^A,g = 0. This implies that to cr^ in eq. 
(|j|) only non-singlet channel partonic subprocesses can contribute. Another 
important property of the coefficient function is that they do not depend on 
the probe 7 or Z so that one can extract the overall pointlike cross section 
or the asymmetry factor A^' . However this is only true if all quark masses 
are equal to zero and if one sums over all quark members in one family 
provided the latter appear in the inclusive state of the partonic subprocess 
(see 0) From (0) we can derive the total hadronic cross section 

^tot(Q 2 ) = W Cdx f 1 dcos9 (x /f \=a T {Q 2 ) + a L {Q 2 ), (10) 
2 ~ Jo J -1 \ dxdcosU J 
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with 



1 r 1 rln H 



(ii) 



where one has summed over all types of outgoing hadrons H . From the 
momentum conservation sum rule given by 



J2 [ dxxDf(x,fx 2 ) = l l = q,q,g 

and eqs. (f|), (|TT|) one can derive 

°k{Q 2 ) = 4ol(Q 2 ) tdxx 
Jo 



(12) 



C Kq (x,Q 2 /^Y + -C Kg (x 1 Q 2 /^ 



(13) 



Finally we also define the transverse, longitudinal and asymmetric fragmen- 
tation functions F^(x, Q 2 )0 



F k H (x,Q 2 



S(Q 2 ) 



a 



dx 



k = (T,L,A). 



(14) 



One observes that the above fragmentation functions are just the timelike 
analogues of the structure functions measured in for deep inelastic electron- 
proton scattering. The calculation of the coefficient functions corrected up 
to order a 2 s proceeds in the same way as done for the Drell-Yan process 
in H and for deep inelastic lepton-hadron scattering in ||. Denoting the 
intermediate vector bosons 7 and Z by the symbol V we have the following 
parton subprocesses. In zeroth order we have the Born reaction 



V 



Y +q 



(15) 



where "/" (/ = q, q,g) denotes the detected parton which subsequently frag- 
ments into the hadron of species H. In next-to- leading order (NLO) one has 

2 Notice that we make a distinction in nomenclature between the fragmentation densities 
and the fragmentation functions F k . 
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to include the one-loop virtual corrections to reaction (|T^) and the parton 
subpro cesses 



V 
V 



V +q + 9 
'g" +q + q. 



(16) 
(17) 



After mass factorization of the collinear divergences which arise in the above 
processes one obtains the coefficient functions which are presented in 
The determination of the order a 2 contributions involves the computation 
of the two-loop corrections to flT5|) and the one-loop corrections to eqs. 
dl6|),(|T7D. Furthermore one has to calculate the following subprocesses 



V 
V 
V 



"q" 



q + g + g 

q + q + g 
q n + q + q + q. 



(18) 
(19) 
(20) 



In reaction (EH)) the two anti-quarks, which are inclusive, can be identical as 
well as non identical. Notice that in the above reactions the detected quark 
can be replaced by the detected anti-quark so that in reaction fl20|) one can 
also distinguish between the final states containing identical quarks and non 
identical quarks. After mass factorization and renormalization for which we 
have chosen the MS-scheme one obtains the order a 2 contributions to the 
coefficient functions which are presented in pj|,[lQ|,|ll||. The most important 
results of our calculations can be summarized as follows. From eq. ( |T3"D and 
the coefficient functions originating from the processes above we can obtain 
o L and <jt corrected up to order a 2 



a L (Q 2 



a. 




2 r 



c n~ y} + CaCf< 



Q 2 24 20231 r T L, Q 2 74' 

~ 11 V - yC(3) + — ] +n f C F T f ^ln- - - 



(21) 
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a T (Q 2 



42 (Q 2 



1 + 



47T 




+ c A c F < 



196 178 
5 SV ; 30 



+ n,C,T / jl6C(3) + - 



(22) 



Addition of cr^ and ax yields the well known answer for cr tot ClTOf) which 
is in agreement with the literature [12]. Hence this quantity provides us 



with a check on our calculation of the longitudinal and transverse coefficient 
functions. Notice that in lowest order cr to t only receives a contribution from 
the transverse cross section whereas the order a s contribution can be only 
attributed to the longitudinal part. In order a 2 both and <jt contribute 
to a tot . 

Because of the high sensitivity of expression ([H]) to the value of a s , the 
longitudinal cross section provides us with an excellent tool to measure the 
running coupling constant. 

To illustrate the sensitivity of <jl on a s we have plotted in Fig. 1 in LO and 
NLO the ratio 



(23) 



vtot(Q 2 ) 

as a function of Q (CM-energy of the e + e~ system). Our input parameters 
for the running coupling constant are = 108 MeV (a^°(M z ) = 0.122) 
and A|| = 227 MeV (af LO (M z ) = 0.118). Fig. 1 reveals that the order a s 
corrections are appreciable and they vary from 48% (Q = 10 GeV) down to 
28% (Q = 200 GeV) with respect to the LO result. Furthermore one observes 
an improvement of the renormalization scale dependence while going from 
to R^ L °. Comparing with the experimental value Rl = 0.057 ± 0.005, 
measured by OPAL 0, one observes a considerable improvement when the 
order a 2 s contributions are included. However there is still a little discrepancy 
between R^ 10 at // = Q = Mz, and the data. This can either be removed 
by choosing a larger Af^g or by including the masses of the heavy quarks c 
and b in the calculation of the coefficient functions. Also a contribution of 
higher twist effects can maybe not neglected (see |J). 

We now want to investigate the effect of the order a 2 contributions to the 
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longitudinal and transverse fragmentation functions Fl(x, Q 2 ) and Fr(x, Q 2 ) 
eq. (|I4|). Here we have summed over the following hadron species i.e. 
H = tt^, K ± ,p,p. Choosing the parametrization of the fragmentation densi- 
ties in [|TJ], with kf = A|| = 190 MeV, we have plotted F[° and Ff LO 
in Fig. 1. We observe that F[° is below the data in particular in the small 
^-region. The agreement with the data becomes better when the order a s 
corrections are included although at very small x F^ LO is still smaller than 
the values given by experiment. This figure reveals the importance of the 
higher order corrections to Fl(x,Q 2 ) for the determination of the fragmen- 
tation densities. We have also shown results for F^ LO and F^ NLO in Fig. 
3 using the same fragmentation density set. Both fragmentation functions 
agree with the data except that F^ NLO gets a little bit worse at very small x. 
Furthermore Fjf and F^ NLO hardly differ from each other which means 
that the order a 2 corrections are small. One also notices that Fl constitutes 
the smallest part of the total fragmentation function F = F? + Fl which can 
be inferred from Figs. 2 and 3. This in particular holds at large x where 
F T » F L . Hence a fit of the fragmentation densities to the data of F T is 
not sufficient to give a precise prediction for Fl and one has to include the 
data of the longitudinal part to provide us with a better set of fragmentation 
densities. This in particular holds for D^(z,fi 2 ) in the small z-region for 
which the order a 2 contribution to Fl will be needed. The latter will cer- 
tainly change the parametrization of the gluon fragmentation density given 
by ALEPH in M and OPAL in §. 
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Figure Captions 



Fig. 1 The ratio R L = cr T /a tot . Dotted lines: R L L °\ solid lines: R% LO . 
Lower curve: /i = 2Q; middle curve: [A — Q; upper curve: \i = Q/2. 
The data point at Q = M z is from OPAL f|. 

Fig. 2 The longitudinal fragmentation function Fl(x,Q 2 ) at fi = Q = Mz- 
Dotted line: Ff °; solid line: Fg L °. The data are from ALEPH § and 
OPAL 0. The fragmentation density set is from fl"3|| . 



Fig. 3 The transverse fragmentation function Ft(x,Q 2 ) at // = Q = Mz- 
Dotted line: F^ LO ; solid line: F^ NL °. The data are from ALEPH @ 
and OPAL |§. The fragmentation density set is from |13 . 
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Figure 1: 
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